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ABSTRACT 

Recent numerical simulations of magnetic reconnection in two dimensions 
have shown that, when the resistivity is strongly localized, the reconnection re- 
gion develops a Petschek-like structure, with the width of the inner diffusion 
region being of the order of the resistivity localization scale. In this paper, we 
combine this fact with a realistic model for locally-enhanced anomalous resistivity 
generated by current-driven microturbulence. The result is a qualitative model of 
the reconnection layer where the size of Petschek's diffusion region and, therefore, 
the final reconnection rate are determined self-consistently in terms of the main 
parameters of the functional dependence of anomalous resistivity on the electric 
current density. We then consider anomalous resistivity due to ion-acoustic tur- 
bulence as a particular case. This enables us to express the reconnection region's 
parameters directly in terms of the basic parameters of the plasma. Finally, we 
apply this reconnection model to solar flares and obtain specific predictions for 
typical reconnection times, which are very consistent with observations. 

Subject headings: MHD — Sun: flares — Sun: magnetic fields 



1. Introduction 



Magnetic reconnection is a basic plasma physics phenomenon of tremendous importance in 
many astrophysical systems (Tsuneta 1996; Kulsrud 1998), as well as in some laboratory 
plasma devices (Yamada et al. 1997), including tokamaks (Kadomtsev 1975; Yamada et al. 
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1994). It has been studied extensively over the past five decades (e.g., Giovannelli 1946; 
Vasyliunas 1975; Biskamp 2000). Historically, the first (and perhaps still the most important) 
application of magnetic reconnection has been to explain the solar flare phenomenon, and it 
is in this context that the earliest reconnection models have been developed. 

The first theoretical model of magnetic reconnection was developed by Sweet (1958) and 
by Parker (1957, 1963). In this model, magnetic field is frozen into the plasma everywhere 
except in a very thin layer where the current density is so high that resistive effects become 
important no matter how small the resistivity is. It is inside this thin current layer that the 
actual breaking (and reconnecting) of the magnetic lines of force takes place, accompanied 
by a violent release of enormous amounts of magnetically stored energy, thus leading to the 
observed flare. The Sweet-Parker theory predicts that the current layer thickness, 5gp, scale 
as <5sp ~ L/y/S, where L is the global system size (typically of order 10 9 cm in the solar 
corona) and S = LVa/t] is the global Lundquist number (here Va is the Alfven velocity and rj 
is the magnetic diffusivity; in the rest of this paper we shall refer to rj as the resistivity). 
Correspondingly, the typical reconnection timescale is found to be r rcc ~ ta(L)\/S, where 
ta(L) = L/Va is the global Alfven transit time. It has been immediately realized that the 
resulting reconnection time turns out to be too long; in the solar corona one typically has 
S ~ 10 12 — 10 14 and ta(L) ~ 1 sec, which leads to r rec of the order of a few months. This is 
in sharp contrast with the typical observed solar flare duration of order 10 2 — 10 3 sec. Thus, 
since its early years, the main thrust of magnetic reconnection research has been to explain 
reconnection rates that are much faster than the Sweet-Parker theory predicts. 

Starting from the 1960s, two major routes toward faster reconnection were proposed. One of 
them was to use the so-called anomalous resistivity instead of the classical Spitzer resistivity 
used in the original Sweet-Parker model (Coppi & Friedland 1971; Smith & Priest 1972; 
Coroniti & Eviatar 1977; Kulsrud 1998). The idea was that, as a current layer forms, 
its thickness becomes so small, and hence the current density becomes so high, that the 
drift velocity of the current- carrying electrons exceeds a certain threshold, such as the ion- 
sound or the electron thermal speed. This leads to the excitation of current-driven kinetic 
microturbulence, which, in turn, provides a more efficient (compared to particle-particle 
collisions) mechanism for the scattering of electrons (via wave-particle interactions). As 
a result, one ends up with a greatly enhanced effective resistivity and, hence, a greatly 
reduced effective Lundquist number. When substituted into the Sweet-Parker scaling, this 
results in a greatly enhanced reconnection rate. In fact, controlled laboratory studies of 
magnetic reconnection have shown a good agreement with a simple Sweet-Parker model 
augmented with some (experimentally measured) anomalously enhanced resistivity (Ji et al. 
1998, 1999). On the other hand, in the solar flare context, the Sweet-Parker model with 
anomalous resistivity gives typical reconnection times of the order of several hours (e.g., 
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Kulsrud 1998), a great improvement over Spitzer resistivity. Among various anomalous 
resistivity mechanisms, the one most frequently quoted has been the ion-acoustic turbulence 
(IAT). Rapid magnetic energy dissipation due to the IAT-driven anomalous resistivity has 
in fact also been invoked to explain coronal heating (Rosner et al. 1978). 

Another possibility leading to shorter reconnection times was proposed by Petschek (1964). 
His very elegant model makes use of a somewhat more complicated reconnection layer geom- 
etry, while still relying on simple resistive magnetohydrodynamics (MHD) without invoking 
any new physics at the microscopic level. The Petschek model actually does not predict a 
unique configuration of the reconnection layer and a unique reconnection rate. Instead, this 
model encompasses an entire one-parametric family of solutions. Each of these configura- 
tions has at its center a small Sweet-Parker-like layer, called the diffusion region, and four 
standing slow-mode shocks emanating from the ends of this central layer. The members 
of this family of solutions can be labeled by the width (or length) A of the inner diffusion 
region. As Petschek noticed, in the Sweet-Parker model the reconnection process had been 
slowed down by the very large aspect ratio L/5 of the reconnection layer. He suggested 
that the width A of the layer's diffusion region does not have to be as large as the global 
size L; if A can be made sufficiently short, the reconnection process will go much faster than 
in the Sweet-Parker model. The maximum value of A corresponds to the Sweet-Parker 
solution, which is, therefore, just one of the family members. The reconnection rate ranges 
from the slowest (Sweet-Parker) rate to the so-called maximum Petschek rate, which scales 
as 1/ log S. We thus see that the relatively strong square- root dependence on the resistivity, 
characteristic for the Sweet-Parker model, is replaced here by a much weaker logarithmic 
dependence; even for a very large S ~ 10 14 , the resulting reconnection timescale turns out 
to be reasonably short, of order 10 2 t a (L). 

This model had remained the favorite model of reconnection until 1980s, when two-dimensional 
(2D) resistive-MHD numerical simulations by Biskamp (1986) showed that, in the case of 
a spatially uniform resistivity, a Petschek-like configuration fails to form and that a long 
(of order L) current layer tends to form instead, consistent with the Sweet-Parker picture. 
This finding has been confirmed in a number of numerical simulations performed by several 
other groups (Scholer 1989; Ugai 1992, 1999; Yokoyama & Shibata 1994; Uzdensky & Kul- 
srud 2000; Erkaev et al. 2000, 2001). A theoretical explanation has been put forward by 
Kulsrud (2001). He noticed that, when the resistivity is uniform, the relatively large trans- 
verse magnetic field that is needed to support the standing shocks in the Petschek model 
is rapidly swept out of the diffusion region by the downstream flow, while its regeneration 
due to the nonuniform merging is not fast enough. As a result, the diffusion region's size A 
(which Kulsrud calls L') increases until it reaches the global scale L and the reconnection 
rate, correspondingly, slows down to the usual Sweet-Parker rate. This explanation has been 
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confirmed numerically by Uzdensky & Kulsrud (2000) (see also Kulsrud 1998). 

It has been noticed, however, that the key assumption leading to the above conclusion was 
the assumption of uniform resistivity (which is a very common assumption in numerical 
simulations in general). Simulations featuring a non-uniform resistivity have shown that 
a Petschek-like structure does form and can be stable whenever the resistivity is locally 
enhanced in some small region near the X-point at the center of the reconnection region 
(Ugai & Tsuda 1977; Sato & Hayashi 1979; Ugai 1986, 1992, 1999; Scholer 1989; Yokoyama 
& Shibata 1994; Erkaev et al. 2000, 2001; Ugai & Kondoh 2001; Biskamp & Schwarz 
2001). A nice plausible theoretical explanation of this phenomenon has again been provided 
by Kulsrud (2001), who analyzed how a locally enhanced resistivity may lead to a more 
efficient regeneration of the transverse magnetic field through non-uniform merging and thus 
to the sustainment of Petschek's shocks. In addition, recent numerical work by Erkaev et al. 
(2000, 2001) and by Biskamp & Schwarz (2001) has shown that the particular Petschek-like 
configuration that forms in the localized-resistivity situation is characterized by the width A 
of the inner diffusion region being of the order of the resistivity localization scale, which 
in this paper we shall call l v (Erkaev et al. 2000, 2001; Biskamp & Schwarz 2001). We 
thus see that resistivity-localization mechanism plays a crucial role in determining the final 
reconnection rate. 

From the point of view of physical reality, the main motivating force behind these localized- 
resistivity studies has been the idea that anomalous resistivity, being such a sensitive function 
of the local current density, may in fact be triggered only in a small neighborhood of the X- 
point, where the current density is highest. In other words, anomalous resistivity can enhance 
reconnection rate not only directly (by simply being higher than the collisional resistivity), 
but also indirectly, via enabling the Petschek mechanism (by being strongly localized). This 
is one of the key ideas behind the so-called spontaneous fast reconnection model suggested 
by Ugai (1986, 1992, 1999; also see Ugai & Tsuda 1977; Yokoyama & Shibata 1994; Ugai & 
Kondoh 2001) on the basis of numerical simulations. 

Up until now, however, with the notable exception of the work by Kulsrud (2001), there have 
been no analytical attempts to combine the Petschek reconnection model with any physically 
realistic model of anomalous resistivity, which would provide unique theoretical predictions 
for the main parameters of the reconnection region, such as the width of the diffusion region 
and the reconnection rate, in terms of the basic parameters of the plasma. The goal of this 
paper is to remedy this situation by attempting to build a simple theoretical framework 
incorporating a particular anomalous resistivity mechanism (ion-acoustic turbulence). 

In § 2 we present our model of the Petschek reconnection layer with a generic form of 
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anomalous resistivity. In particular, in § 2.1 we describe the basic elements of the model, e.g., 
the Sweet-Parker relationships for the inner diffusion region and the functional dependence 
of anomalous resistivity, r], on the current density j. In § 2.2 we discuss possible solutions 
and explore their stability. In § 2.3 we calculate the reconnection rate in terms of the model 
parameters. § 3 is devoted to the specific case of anomalous resistivity due to the IAT; in this 
section we make use of the well-developed theory of ion-acoustic turbulence to express all 
major reconnection layer parameters, including the reconnection rate, in terms of the basic 
plasma parameters (such as the magnetic field strength, plasma density, and the electron 
and ion temperatures). We apply the obtained results to solar flare environment in § 4 and 
find a very reasonable agreement in terms of general timescales. We list our conclusions 
in § 5. 



2. Model of the Reconnection Layer 

2.1. Three Main Ingredients of the Model 

We now describe our model of a Petschek-like reconnection configuration that is formed in 
the presence of anomalous resistivity due to a current-driven microturbulence. This model 
is semi-empirical and does not pretend to describe any real physical system in full detail. 
We believe, however, that it correctly captures the most critical qualitative features of the 
reconnection phenomenon. 

Our model represents a synthesis of the following three external ingredients that are com- 
bined to build a complete and, hopefully, self-consistent description of the system: 

1) the numerically-observed fact that, whenever the resistivity is strongly localized, a Petschek- 
like structure tends to develop, with a characteristic width of the central diffusion region 
being of the order of the resistivity localization scale (Erkaev et al. 2000, 2001; Biskmap & 
Schwarz 2001); 

2) Sweet-Parker model for the central diffusion region of the Petschek configuration; 

3) a physically motivated model for the anomalous resistivity. 

We shall now proceed with the description and integration of these three main ingredients. 
Let y be the direction along the layer and x across the layer. Figure 1 shows schematically 
the central part of the Petschek-like configuration discussed in this paper. The shaded 
rectangular area (of characteristic thickness 5 and width A) in the center of the Figure is the 
inner diffusion region where the electric current density is concentrated. This is the region 
where the effective plasma resistivity plays an important role and where the actual breaking 
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of magnetic field lines takes place. 

For definiteness, we use Ampere's law to define the thickness 5 of this central current layer 
in terms of the central current density j = j(0, 0) and the outside magnetic field B as 

(1) 

In addition, suppose that along the midplane x = the current density has a certain pro- 
file j(y)- We define the characteristic width A of the diffusion region in terms of the func- 
tion j(y) as the distance from the center y — to the point where the current density drops 
by a factor of e: 

j(y = A)= Jo /e. (2) 

Similarly we define the characteristic scale l v for the variation of the plasma resistivity r](y) 
along the layer as 

v(y = Q = v(y = o)/e. (3) 

In general, the two scales A and l v may be very different. However, as discussed in the 
introduction, recent numerical simulations have shown that, if the resistivity is strongly 
localized (l v <C L, where L is the global system size), a Petschek-like configuration develops 
with the width of the inner diffusion region being of order of the resistivity localization scale: 
A ~ Ijj. This important condition serves as the criterion for selecting a unique solution out of 
the entire family of Petschek configurations. We shall call this unique Petschek configuration 
an equilibrium configuration. In our present analysis we shall rely on this empirical finding 
as on one of the important building blocks of our model. 

Note that, in general, A and l v in this equilibrium Petschek configuration may differ by a 
finite factor, so let us define a dimensionless parameter K such that the equilibrium configu- 
ration has A = A eq = Kl v . In our analysis it will actually be more convenient to use another 
dimensionless parameter to describe the equilibrium Petschek configuration. We shall call 
this parameter £ and define it in terms of j(y) as jil^) = £jo/ e - The parameters K and £ 
are related to each other, but the exact relationship depends on the detailed profiles of 7]{y) 
and j(y); it is, however, unimportant for the purposes of the present paper: we shall only use 
the fact that, since j(y) is a monotonically decreasing function of y, we have £ > 1 whenever 
A cq > l v (K > 1) and £ < 1 whenever A eq < l v (K < 1). We expect £ to be of order one, its 
precise value depends on the details of the problem and must be determined from numerical 
simulations. The analysis in this paper will be restricted to a situation where A eq , i.e., the 
equilibrium value of A, is somewhat greater than l v ] thus we take £ > 1. This choice is made 
purely for reasons of convenience, as will be elucidated below. 
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We shall now describe the other two main components of the model. The second ingredient 
is the model for the inner diffusion region. In Petschek's theory, this region is a Sweet- 
Parker-like current layer with the thickness 5 related to its width A via the relationship: 

{ = (4) 

where 

S A , — (5) 

Vefi 

is the Lundquist number for the scale A. Here we shall take r] e s = r](j ), the resistivity at 
the very center of the layer. Then we can rewrite the above expression for S as 



This is a very important relationship that we are going to invoke many times throughout 
the paper. 

The next question is what actually determines the value of and hence A. We suggest 
that these values have to be determined by the properties of the function r](j). We thus 
need to introduce the third ingredient, namely, a physically plausible model for anomalous 
resistivity, expressed in terms of a single function r](j). 

Our choice of the function 77(^7 ) is motivated by anomalous resistivity due to a current- driven 
microturbulence. Correspondingly, we here adopt a very simple, minimal model for 
which, however, has to exhibit the following general properties: first, there exists a current- 
density threshold, j c , for triggering anomalous resistivity, and second, the rapid rise of i] after 
the threshold is exceeded stops at some large but finite value rji, after which continues 
to rise with increased j more slowly. 

Thus we take r\ (j) to be a prescribed function whose behavior can be described as follows 
(see Fig. 2): 

1) For j < j c , r](j) is constant and equal to a small collisional resistivity i] . 

2) At the critical current density, j = j c , 77 rises rapidly, essentially jumps to a value rji ^> r] . 
For our convenience, we actually introduce some small finite width to this jump, Aj <C j c , 
and assume that r](j) is a linear function in the interval j G [j c ,jc + Aj]. The exact value of 
Aj is unimportant in our analysis. 

3) When j is increased even further, j > j c + Aj, the resistivity continues to rise monoton- 
ically with increased j, but in a much slower manner. 1 In particular, we take r\ ~ j in this 



1 Numerical simulations sometimes assume a stronger saturation of anomalous resistivity by imposing a 
strict upper limit on 77 (Yokoyama & Shibata 1994). 
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region, which is the case for ion- acoustic turbulence (e.g., Bychenkov et al. 1988). Thus we 
take 

v (j> j c + Aj)^ Vl ^. (7) 

Jc 

We see that the dependence of the anomalous resistivity r\ on the current density j involves 
three parameters, j c , t]q, and rji. In addition to these parameters, it is convenient to introduce 
the critical layer thickness 5 C as a derived quantity: 

4ttj c 



For the analysis in this section this generic level of description will suffice. In § 3, however, we 
shall consider a very important specific example of anomalous resistivity due to ion-acoustic 
turbulence and will give specific expressions for these parameters. 

Here, however, we wish to make a remark concerning the manner in which we are going 
to use the above prescription for r](j). In our model, rj is determined solely by the local 
value of j. We realize, of course, that in reality the coefficients j c , r]o, and rji will depend on 
the plasma density and temperature, and so will vary from one point to another inside the 
reconnection layer. In our simple reconnection model we shall, however, ignore this aspect 
and assume that these parameters are constant in space and time and hence the same profile 
v(j) applies everywhere. In other words, r/ is determined solely by the local value of j. 



2.2. Two possible solutions and their stability 

We shall now use the 7](j) dependence shown in Figure 2 to find the correct value of jo 
that corresponds to the equilibrium Petschek solution. First, the requirement that A = Kl v 
can be interpreted as follows. At y — 0, let us have some value jo an d the corresponding 
value i](j ). The requirement that % = j{l v ) = £jo/ e means that r](y = l v ) = r/^jo/e). But, 
by definition of l v , i](y = l v ) = r](j )/e. Thus, jo is determined using function i](j) from the 
condition: 

vUo) = evU*) = eviUo/e) ■ (9) 

The solution of this equation for a given £ can be found by drawing a set of parabolae rj oc j a 
where a — 1/(1 — ln£), and selecting out of this set those parabolae for which the points 
of their intersection with the curve r](j) are separated by a factor of e in their values of 77, 
as shown in Figure 3. It is clearly seen from this Figure that, for any £ > 1 (and, hence, 
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a > 1) and for the general shape of r)(j) adopted in this paper, one finds, in principle, two 
such solutions. 

The first solution (curve I) corresponds to jo lying on the rapidly rising part of the rj(j) 
curve, il e [j c Jc + Aj], while j{ = ^{Q = £j /e < j c . Then, r] T (y = l v ) = r] and 
hence ?7 J (jo) — e Vo- We thus see that the corresponding resistivity enhancement in the inner 
diffusion region is not very large in this case, just by a factor e over the collisional value 770- 

The second solution (curve II) boasts much higher resistivity enhancement. In this solution, 
jl 1 falls within the narrow rapidly rising part of the 77 (j) curve, and hence 

3" = VWJ ^ ejc/t ■ (10) 

Correspondingly, the resistivity at the center of the diffusion region is 

v"tio)=vti")=rne/Z^r i0 - (H) 

Any Petschek-like configuration with j between jl and j^ 1 , and hence with erjo < rj < er)i/£, 
will not be in a steady equilibrium and instead will evolve so as to increase j . Indeed, if 
vUo) > er 1o, then 77 ( ) > i] and hence > j c . But then, since j < ej c /£, we have 
j*/jo > e /C an d therefore A > Kl v . As we know from numerical simulations, this leads to 
shrinkage of A. As we demonstrate below, this shrinkage, in turn, results in a decrease in 5 
and hence in an increase in jo- The evolution will then presumably reach a stationary state 
when (the stable) solution II is reached. 

Now, which one of the two solutions I and II will be realized in a real physical system? We 
suggest that, in general, the answer to this question is determined by the stability of the 
solutions with respect to a small change in jo. In particular, we shall demonstrate that the 
first solution is unstable while the second one is stable. This will enable us to conclude that 
the system will evolve towards the second solution corresponding to higher rj and hence to 
a higher reconnection rate. 

We shall first demonstrate that the stability properties are largely determined by the relative 
size of the slopes of the function r](j) at j = jo and j = j* = £jo/ e - We start by presenting a 
rather general stability analysis that works for any monotonically increasing function 
Hence, in order to apply this analysis to our solution I, we need to modify our r)(j) slightly 
by assuming that this function has a non-zero positive (although arbitrarily small) slope 
v'U) > in the region j < j c . We then discuss a somewhat modified treatment for the 
case when rj' (j) is exactly zero in this region. These modifications are not essential and the 
conclusion regarding solution I being unstable is the same. 
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So, first let us assume that r/(j) ^ but has a small positive value at j < j c . Let us 
consider an equilibrium Petschek-like configuration with A = Kl rj . Let this configuration be 
characterized by unperturbed values 5, j , j*, A, l v , etc. Now, imagine that at t — this 
equilibrium Petschek configuration is suddenly perturbed while preserving its Petschek-like 
character, that is changed into a neighboring Petschek configuration. This new Petschek 
configuration is, generally speaking, not in equilibrium, i.e., it does not satisfy condition (9). 
It will then evolve through a sequence of Petschek states. For simplicity, here we envision the 
process as occurring in two parts with different timescales: the adjustment between A and S 
to conform to the Sweet-Parker structure of the diffusion region is instantaneous, whereas 
the adjustment of A to l v occurs on a longer timescale. In reality this might not be true and 
these two processes may occur on the same time scale (namely, the Alfven crossing time). 
But here we are only interested in the direction of the evolution, i.e., whether the perturbed 
system will move away or towards the equilibrium, so our qualitative analysis should be 
sufficient. 

In particular, let us imagine that at t — the thickness S of the diffusion region is reduced 
slightly and hence the central current density is correspondingly slightly increased: 

5 -> 6 = 6(l-e), (12) 
jo -> jo=jo(l + e), (13) 

where e < 1. 

An increase in j by a factor (1 + e) leads to an increase in the resistivity at the center: 

V(io) -> vUo) = V(jo) + v'Uo)joe ■ (14) 



For the system to remain a valid Petschek solution, we require that the inner diffusion region 
remains a Sweet-Parker layer at all times, so the width A of the current distribution will 
automatically adjust to a new value, A(0), which, according to equation (6), is 



A -> A(0) = V A 



vUo) 



e(2 + 



Vtio) 



(15) 



This equation expresses the initial reduction in A in direct response to the reduction in 5 
and hence to the related increase of jo and r)(jo) in the Sweet-Parker model for the diffusion 
region. For example, if we are considering solution I, j e [j c , j c +Aj], then r)'(j ) ~ rji/Aj ^> 
Vi/jo > v(jo)/jo, an d hence the last term is dominant, i.e., the initial change in A is mainly 
due to the change in resistivity: 



A'(0)*A(l-c^^W(l 
v VVo) f 



vUo) Aj' 



(16) 
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Similarly, if we are considering solution II, then j > j c + Aj, and r\ oc j so that 7/ (jo) = 
vUo)jo- Then 

A n (0) ~ A^l -3e) . (17) 



Let us now ask what the new value l v of the width of the enhanced-resistivity region is and 
how it compares with the current sheet width A(0). 

From equation (14) we see that the point y = l v moves to a location where r](y = l v ) is 
increased by rj'(jo)joe/e over its unperturbed value. On the other hand, this increment is 
equal to ?/(j*)(j* ~ 3*)- Therefore, j* has to increase by the amount 



3* 3* 



3o (jo) = £ (jo) 
e rj' (j*) f 77' (j*) 



(18) 



We have, by combining the expressions for the changes in j and j*, 



M-) [ i+e (m-0 



(19) 



We are now going to use this result to evaluate the change in L 



Because of the symmetry with respect to y — 0, the ratio j(y)/jo should be an even function 
of y with the characteristic scale A: 



j(y) 

jo 



= F(Y) , 



where 



y = 4-. 



y_ 

A 



(20) 



(21) 



Let us assume for simplicity that during the initial perturbation the shape of this function 
does not change. We can then write (here = l v /A) 



t = F(Y t ) = F(Y t ) l-nfr-Y, 
Jo 



where 



^ F> { Y) 



F{Y) 



= -|F'(y,)>o. 

Y=Y, £ 



(22) 



(23) 



By comparing this with the above result (19) we see that the change in Y* is equal to 



' V(jo) 



(24) 
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This equation tells us whether the ratio Y* = /^/A increases or decreases from its equilibrium 
value 1/K. We see that the result depends on the relative sizes of the slopes of function r)(j) 
a t j = jo an d j = j*. In particular, in the case of solution II, r) f (j ) < £?/C?*)> an d then it 
follows from equation (24) that l v is reduced by a lesser factor than A. As a consequence, 
the resulting perturbed Petschek configuration has A < A cq = Kl v . We then expect, based 
on the results of numerical simulation mentioned above, that in this case A will increase 
towards A eq , thus negating the effect of the initial-perturbation decrease in A. 

In contrast, if i]'(jo) > Cv'U*) > 0> then l v is reduced by a larger factor than A; thus 
the resulting perturbed Petschek configuration has A > A cq = Kl v . Invoking again the 
results of numerical simulations, we conclude that the system will evolve in such a way as to 
decrease A further. In fact, this scenario describes what happens in the case of solution I, 
but its application to this solution involves some subtlety. Indeed, notice that one can use 
equations (19) — (24) only if ?/(j*) ^ 0. Thus, in order to apply this general analysis to 
solution I, we need to modify the function slightly to make it monotonically increasing 
in the region j < j c . In other words we tentatively assume that this function has a non-zero 
positive (although arbitrarily small) slope > in this region. This modification is not 
essential; it needs to be introduced here purely for our technical convenience, to be able 
to use our general equations (19) — (24). We do not actually need to rely on it in order 
to derive the same conclusions regarding solution I. Indeed, if we insist on our unmodified 
function r](j) with the slope r)'(j) being exactly equal to zero for j < j c , then, instead 
of equations (19) — (24), we can use an alternative argument to describe the response of 
solution I to the initial perturbation. This argument is in fact even simpler and physically 
more transparent than the one described above; it goes as follows. 

Since jo in solution I lies on the steeply rising part of the i](j) dependence, the central 
resistivity is increased by a relatively large amount given by equation (14). On this part of 
the curve, we have r]'(j ) ~ r/i/Aj, and so 

V(jo) = er] + erj^ . (25) 
Aj 

Then, r](l v ) = r}(j )/e = rj + ^(jJAj), and, therefore, j m G \j c ,j ], i.e., j* m j . If 
the perturbed current density j(y) has a characteristic scale A, then j* pa j implies that 
I, <A. Thus we see that, as the principal effect of the initial perturbation, the size l v of the 
resistivity enhancement region drops very sharply, whereas all other quantities change rather 
smoothly. Then we can again invoke our empirical fact that a configuration like this will 
not stay stationary but will evolve so that the current layer width A will tend to decrease 
to become comparable with l v , in agreement with our expectations. 

But as A changes during the subsequent evolution for t > 0, there will be a feedback on 8. 
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Indeed, since the diffusion region always remains a Sweet-Parker layer, 5 and A are always 
related via equation (6). Since r](j ) is a monotonically increasing function of j oc 1/5, 
we see that 5 is a monotonically increasing function of A. For example, for solution I, 
we can, using the fact that jo is always very close to j c , approximate 5 by 5 C everywhere 
except where it appears in the combination (5 C — 5). We thus can express the resistivity as 
7](j ) ~ i]i(j c / Aj)(5 c — S)/5 C and, substituting this expression into equation (6), obtain 

( Sc-h 1 _ 5 c V A 5 c Aj 

v s c J A m Jc ' 

We thus see that the decrease in A(t) leads to a proportional increase in 5 C — 5(t). In a 
similar manner, for solution II we find, using equations (6) and (7), that 5 oc A 1 / 3 . In either 
case we see that after the initial perturbation, the subsequent evolution of the system is such 
that 5 changes in the same direction as A. 

Now we are ready to apply the above general results to the issue of stability of our solutions I 
and II. We see that in the case of solution I, the initial decrease in 5 and A tends to amplify 
further. This means that this solution is unstable: a slight increase in jo will lead to further 
increase and thus the system will move away from the initial equilibrium, towards higher 
values of central current density and resistivity, until it approaches solution II. 2 Solution II, 
on the other hand, is stable; indeed, as have shown above, after its initial decrease [equa- 
tion (17)], A and hence 5 will tend to increase, thereby reducing the initial perturbation. 
Thus, in this case the feedback is negative and hence solution II is stable. 

It is important to realize that, in practice, solution I is likely to be irrelevant for another 
reason (in addition to being unstable). Indeed, real astrophysical systems, including solar 
flares, always have a finite (albeit very large) global size L. For the Petschek model to 
work, we must have A < L. This condition imposes an upper limit on A and hence, via 
equation (6), a lower limit on rj(jo). For example, for j on the rapidly rising part of the rj(j) 
curve, this condition can be cast as 

VUo) > Vmin(L) = Si . (27) 

If ^min > e Vo, then solution I does not even exist for a given global size L. For typical solar 
corona conditions, collisional resistivity i] is so small that inequality (27) is not satisfied; in 
other words, one gets a very large value for A 7 , much larger than L ~ 10 9 cm. 



2 At the same time, a slight initial decrease in jo will presumably lead to a further decrease in jo and 
hence an increase in A, until the system reaches the stable Sweet-Parker configuration with A = L and 
77 ( jo ) = Vo- This scenario, however, can only occur if the system's size L is tremendously large; it is therefore 
of no importance to solar flares, as discussed below [see equation (27)]. 
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Here is how the requirement that r)(jo) > Vmm(L) comes into play from the evolutionary 
point of view. As a reconnection layer starts to form, its thickness decreases rapidly until it 
reaches 5 C ; after that, 5 and the central current density jo both stay approximately constant 
while 77 ( jo ) increases rapidly. When 77 ( j ) reaches ?7 m i n , the corresponding width A of the 
diffusion region, determined by equation (6), becomes equal to L. From this moment on, 
the reconnection system can be described by the Petschek model. The particular Petschek 
configuration found exactly at this moment (i.e., when 77 = 77 min and hence A = L) is that 
of Sweet-Parker. If r? min > er? , then j(l v ) > j c and hence is very close to jo! this means 
that resistivity is very strongly localized, l v <C A. Therefore, this configuration will not be 
at equilibrium, and the system will evolve through a sequence of Petschek configurations 
with ever-increasing jo and Each of these configurations with jo on the rapidly rising 

part of the 77 (j) curve will have A > A cq = Kl v until the configuration corresponding to the 
stable equilibrium solution II is reached. 

We have thus demonstrated that a physical system with anomalous resistivity of the type 
shown on Figure 1 will evolve towards a stable solution II, which is characterized by the 
following values of the central current density and resistivity: 



The exact coefficients in these expressions require a precise model for r)(j), etc... 

Finally, let us make a couple of remarks regarding our choice of £. In the above considerations 
we have assumed that £ > 1. Notice that the case £ = 1 (A eq = l v ) is, in a certain 
sense, degenerate: in addition to an unstable solution I, one gets a continuum of neutrally 
stable solutions II that correspond to j > ej c . Our present analysis does not allow us to 
discriminate among these solutions; they all appear to be equally plausible. If, however, one 
considers the unique stable solutions II with £ > 1 and takes the limit £ — > 1, then one 
arrives at a single limiting solution (having jo = ej c ). Thus, we use this limiting solution to 
extend our family of solutions with £ > 1 to include the £ = 1 case. 

However, if £ < 1 (A cq < l v ), then our present analysis does not work: one still has an 
unstable solution I, but no solutions of type II exist, at least as long as 7?(j) remains a linear 
function. This suggest that the system will continuously evolve towards higher and higher 
values of jo and 7?(j ) and hence the reconnection process will continuously accelerate until 
some new resistivity saturation mechanism sets in. In this case, the system will stabilize and 
a steady state will be achieved at a much higher level of jo than the level jo ~ j c discussed 
here. 



Jo ^ ejc/Z, 



(28) 



(29) 
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2.3. Reconnection Rate 



Thus, we now have a stable steady state Petschek-like configuration with the central diffusion 
region characterized by the thickness 5 ~ 5 c /( and the central resistivity rj = (rji. Here, ( 
is a finite constant; according to our analysis, ( = e/£. However, we realize that the simple 
qualitative model presented here is not adequate for providing any numerical accuracy; it 
should only be used for order-of-magnitude estimates. Therefore, in the following we shall 
just keep using the above expressions for 5 and rj(jo), parameterizing our ignorance by a 
constant ( = 0(1). 

Let us now ask, what is the reconnection rate associated with this configuration? Using 
expression (4) we find the aspect ratio of the central diffusion region: 

A ™ S - (30) 



where we define 

S.= ™. (31) 

m 

[Note: the parameter S 1 * defined here is completely different from Kulsrud's (2001) S*.] 

The reconnection velocity is (assuming that the density stays roughly constant inside the 
reconnection layer): 

^4 = C 2 S,-, (32) 
and the typical reconnection timescale is 

L 



t a (L)SJC, (33) 



V 

v rec 

where r A (L) = L/V A is the global Alfven crossing time. 

Note that the aspect ratio (30) and hence the reconnection velocity (32) turn out to be 
independent of the global system size L, but only depend on the local plasma parameters. 

Also note that our expression (32) for the reconnection velocity differs from the result 
V iec /VA ~ (8 c r]* /VaL 2 ) 1 / 3 obtained previously by Kulsrud [see equation (26) of Kulsrud 
2001]. We believe that this discrepancy can be attributed, at least partly, to a somewhat 
different functional form of the ?7(j)-dependence adopted in his paper. 
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3. Petschek Reconnection in the Presence of Anomalous Resistivity due to 

Ion-Acoustic Turbulence (IAT) 

In this section we assume that the anomalous resistivity enhancement is caused by the 
scattering of the current- carrying electrons off ion-acoustic waves, which are themselves 
excited by the ion-acoustic instability when the current density exceeds a certain threshold. 
The theory of ion-acoustic turbulence and the associated with it anomalous resistivity has 
greatly progressed over the past 40 years (e.g., Kadomtsev 1965; Rudakov & Korablev 1966; 
Sagdeev 1967; Tsytovich & Kaplan 1971; Biskamp & Chodura 1972; Coroniti & Eviatar 
1977; Bychenkov et al. 1988). This theory is now very mature and seems to be capable of 
producing reliable quantitative results regarding anomalous resistivity. In this paper, we use 
the results presented by Bychenkov et al. (1988). 3 

Following the analysis of ion-acoustic turbulence by Bychenkov et al. (1988), we adopt 

j c = a x en e v s , (34) 

where v s = ^/ ZT e /rrii is the ion sound speed, and ai = 0(1). In particular, according to 
equation (2.143) of Bychenkov et al. (1988), a x ~ 2.14. 

Further, we compute anomalous resistivity > j c + Aj) by using the expression (2.148) 
of Bychenkov et al. (1988) for a = o- anom (E) to express i] = c 2 /4ira in terms of the electric 
current density j = aE. We thus obtain 4 



VU) ~ 4tt-0.16 u 2 e \\ m ) V8™ e T e ~ 36 u pe % V m en e v s ' 1 ) 

We then estimate rji by extrapolating this dependence down to j = j c + Aj ~ j c . Using 
equation (34), we get 

/ ■ \ ° 2 ZT e I Zm e 

Vi = VUc) =a 2 —\ , (36) 

Up e Ti V rn,i 

where 

a 2 ~ H ~ 0.06 . (37) 



3 We note, however, that here we use the results of the theory of ion-acoustic anomalous resistivity that 
has been developed for a homogeneous plasma without magnetic field. We acknowledge that the resistivity 
may be modified by both the presence of the magnetic field and by the fact that in our analysis turbulence 
is presumed to be confined to the very small central diffusion region, and is thus strongly inhomogeneous. 

4 Note that the value for anomalous conductivity given by Bychenkov et al. (1988) differs by a numerical 
factor of order one from the famous Sagdeev's (1967) formula. 
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For a pure hydrogen plasma (Z — 1, mi — m p ), we have: 



^ = 1.4.10- s — ^. (38) 

As for 770, we take it to be the classical collisional Spitzer resistivity 

c 2 0.02A c 2 , . 

Vo = Vs P = ~ A , 39 

47T(TSp N D U pe 

where A is the Coulomb logarithm and 

N D = n e X 3 De , (40) 



Xoe = \fT e /Ann e e 2 being the electron Debye radius. We see that we may expect a potential 
resistivity enhancement on the order of 

r?i N D T e lm e , . „ , . 

— = -ZTT\— »!■ 41 

2A Tj V mi 

For this model of anomalous resistivity, 5 C and 771 are given by expressions (8) and (36) and 
then we can express 5* as 

01^2 cf s Zi e Zm e 

Expressing the ratio Va/^s in terms of the composite electron plasma beta parameter, 

A-^, (43) 



as Va/^s = y/2/~/3l, we can rewrite (42) as 

c YA Ti m l AA\ 

cvPe ZL e Zm e 

where 

Pi 

a 3 = — ~ 11 , (45) 
aia2 

where we substituted a\ = 2.14 and made use of equation (37). 
Note that from equations (8) and (34) it easily follows that 

5 r 1 [2 



CLi di V Pe 
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where 

di = — (47) 

is the ion skin depth. Substituting this very useful expression, along with equation (44), into 
our equation (30) for A, we find 

A ^e d 'z^zr e ^- < 48 > 

where 

a 4 = -5— ~ 7.3 . (49) 
afa 2 

We also get an expression for the Alfven crossing time for the diffusion region — one of the 
most important timescales in the problem: 

^ A ) = ^C^& (50) 

In addition, substituting equation (44) into equation (33) for the reconnection time r rcc , we 
find a very simple relationship expressing r rec in terms of the light crossing time L/c and j3 e : 

a 3 L Ti rrii 

(51) 



C 2 ^W e c ZT e Zm e 



For a pure hydrogen plasma (Z — 1, mi — m p = 1836 m e ), we get 



, 4 V A Ti 



S* ~ 2-10 4 — ^=-f, (52) 



- e 



A ~ i. 3 .10 4 r^^-f , (53) 



r A (A) ~ I.S-IOV^-V^ 1 , (54) 
C 2 v^ c T e v ; 



Note that in all these expressions n e and T e are to be taken at the center of the reconnection 
layer (x — y — 0), while the magnetic field _B is the reconnecting magnetic field outside the 
layer, at x > 5, y — 0. 
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4. Application to Solar Flares 

Let us now try to apply the above results to typical solar flare conditions and see whether 
our model is able to explain the very short time scale of impulsive flares. In order to be able 
to make quantitative estimates, we shall first need to discuss the values of some relevant 
plasma parameters. 

Table 1 lists the values of the key parameters of our model, along with the values of some 
fundamental plasma parameters, for two sets of conditions. Both sets are calculated for fully- 
ionized pure hydrogen plasma. The first set (column III) illustrates the fiducial solar coronal 
conditions: B = 100 G, n e = 10 9 cm" 3 , and T e = T { = 2 • 10 6 K ~ 200 eV. The second set 
(column IV) corresponds to the fiducial solar flare conditions which will be discussed below. 

For the parameters in column III we see that the characteristic reconnection time turns 
out to be no more than an order of magnitude longer than the observed flare duration 
time (which is typically of order 10 3 sec), but it is, apparently, still not sufficiently short. 
Notice, however, that the fiducial solar corona parameters used in column III may not be 
appropriate for the center of the solar flare reconnection layer. Indeed, one can expect that 
the turbulence will lead to rapid heating of the plasma, resulting in a substantial rise of the 
electron temperature. 

We can then ask whether the electron thermal pressure at the center of the current layer will 
grow to a level where it becomes comparable to the outside magnetic pressure, f3 e = 0(1). 
This is a very important question because, according to equation (46), the value of f3 e controls 
the regime that the system finds itself in. In particular, if 5 < di, then the diffusion region 
needs to be described in terms of electron MHD (or Hall MHD), the theory of which in the 
reconnection context has recently been greatly advanced by a number of researchers (e.g., 
Drake et al. 1994; Biskamp 1997; Bhattacharjee et al. 2001). Thus, we see that (3 e is a very 
important parameter, whose value may have a profound influence on the applicability of the 
anomalous resistivity model adopted in this paper. Let us now ask what value for f3 e one 
can expect in the case of magnetic reconnection in solar corona. 

First, notice that it seems inevitable that, if there is no axial (or guide) magnetic field 
component, that is when B z = 0, then one has to have f3 e = 0(1). Indeed, in this case 
the basic requirement of pressure balance across the current layer dictates that the plasma 
pressure at the center of the layer's diffusion region be equal to the pressure of the magnetic 
field outside of the layer. This means that (5 e — 1 when T e ^> Tj (and hence the ion pressure 
inside the layer is negligible) and j3 e = 1/2 when Tj = T e (and hence the ion pressure is 
equal to the electron pressure for pure hydrogen plasma). This conclusion holds regardless 
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of the energy budget balance. For example, if there is no effective cooling mechanism, then 
the electron density n e does not change significantly (i.e., by more than a factor of order 
one), whereas the temperature increases up to the "equipartition level", T eq = Bl/%Tcn e , 
which is about 3 • 10 s K for n e = 10 10 cm" 3 and B = 100 G. On the other hand, if there 
is some effective cooling, the plasma temperature cannot reach such a high value and then 
the density is increased instead to maintain the pressure balance across the layer. In either 
case, one finds (5 e = 0(1) and hence, according to equation (46), 5 ~ d^. This means that 
the system will require a Hall-MHD description as soon as, or even before, the anomalous 
resistivity becomes important. 

In the case of solar flares, however, it is unclear whether this situation is present. Indeed, 
there exist a possibility for maintaining the pressure balance across the layer with /J e < 1. 
This scenario requires two things, both likely to be relevant in solar flare environment: the 
presence of a non-zero axial component of the magnetic field B z and some plasma cooling 
mechanism. Indeed, when the plasma is cooled efficiently, the increase in thermal pressure 
that is required to maintain the pressure balance cannot come from the increase in tempera- 
ture, and so the plasma tends to compress inside the reconnection layer. This, in turn, leads 
to the proportional compression of the guide field component (here we are neglecting the 
resistive decoupling between the guide field and the plasma): 

B z inside inside /rn\ 

FM = 7T\ • (56) 

- 2 | outside ' | outside 

(Here the subscript "inside" corresponds to the center of the reconnection layer, x = y = 0, 
whereas the subscript "outside" corresponds to the plasma above the layer, y = 0, x ^> 5.) 

If cooling is so strong that the resulting central temperature is small compared with the 
equipartition temperature, then the plasma pressure can be neglected in the pressure balance. 
The pressure balance is then achieved with the increased guide field pressure inside the layer 
balancing the reconnecting field's pressure outside the layer: 

-B 2 |insidc B z |outsidc = Bq . (57) 

Thus, -B z |i n side is determined from the pressure balance (57) and then the ratio B z \i ns [^ c /B z | outS idc 
determines the compression factor and hence the central density n e |i ns id e - Typically one might 
expect -B 2 |outside ~ Bq, and so B z (and hence n e ) is increased at the center of the layer by a 
factor of order one. As for the central electron temperature, it is going to be determined by 
the balance between the turbulent ohmic heating and the cooling due to electron thermal 
conduction. 5 This is an important and very complicated problem and its detailed treatment 



5 In the context of solar flares the radiative cooling of the diffusion region, including both the 
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lies outside the scope of this paper. Therefore, here we shall give only some very simple 
estimates. 

The characteristic ohmic heating time can be evaluated as 

Theat ~ —q- , (58) 

where the ohmic heating rate per unit volume roughly is 

Jo 2 BZ V (j ) Bl 1 

Q ~7~4^^~4^(A)' (59) 

so that 

r hcat ~r A (A)(3 e , (60) 

where we have dropped numerical factors of order unity. We thus see again that, in the case 
where there is no cooling, a cold fluid element entering the reconnection layer will be heated 
up very rapidly; in fact, the total time that the element spends (and is being heated) inside 
the inner diffusion region [of order 7a (A)] is long enough for the plasma thermal pressure to 
reach (5 e = 0(1). 

If there is efficient cooling due to electron thermal transport, then an absolute lower bound 
on the electron cooling time is set by the time it takes a freely streaming thermal electron 
to leave the inner diffusion region of size A: 

T ~" ~ it ~ Ta[A) ^wM ■ (61) 

^th,e VPe V m i 

Then, equating r c ™° and 7h ea t, we get a lower bound on f3 e : 



This is the regime that is illustrated in column IV of Table 1. Here, we use the following values 
of plasma parameters for our "fiducial solar flare conditions": B = 100 G, n e = 10 10 cm -3 , 
T e = 3 • 10 7 K ~ 3000 eV, and T { = 3 • 10 6 K ~ 300 eV, which correspond to f3 e ~ 0.1. 
The resulting reconnection time scale r Tec is of the order of a hundred seconds, which is fast 
enough to explain the observed very short duration of the impulsive phase of solar flares. 



bremsstrahlung and cyclotron mechanisms, appears to be ineffective, as the characteristic radiative cool- 
ing time is much longer than the time Tyi(A) that a fluid element spends inside the diffusion region. 
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An upper limit on r coo i is given simply by the time that a fluid element spends inside the 
inner diffusion region, i.e., the Alfven transit time ta(A). This corresponds to the upper 
limit /3™ ax = 0(1), just as in a situation without a guide magnetic field or without cooling. 
The electron temperature then grows to about the equipartition value T cq ~ 3 • 10 8 K (for 
n e = 10 10 cm' 3 and B = 100 G). 

Here are a few more numbers: 
23, 45.32, 18650, -0.652. :) 

Thus, we can constrain f3 e to lie between (3™ m ~ 10 _1 and f3™ ax = 0(1). We then see from 
equation (46) that even in the case of the lowest possible (3 e , the thickness 5 of the inner 
reconnection layer is roughly of the same order of magnitude as di. This suggests that the 
Hall-MHD regime is likely to be at least marginally important in the physics of solar flares. 
In this case the nature of anomalous dissipative processes may differ from the simple effective 
resistivity rj(j) due to ion-acoustic turbulence as described in § 3 and thus a more elaborate 
theory is needed. 

Finally let us address one more question related to the applicability of the ion-acoustic 
regime. One can raise the objection that, in order to excite ion-acoustic instability, the 
condition T e ^> T« needs to be satisfied; if instead the plasma is nearly isothermal, with 
T e ~ Tj, then Buneman instability (Buneman 1959) can in principle be excited, but at a 
higher current-density threshold, j® micman ~ en e w t h,c ^ en e t> s . To address this problem, let 
us consider a plasma that initially (i.e., before the onset of reconnection) is isothermal. As the 
reconnection current layer is been formed, the current density in the layer gradually increases 
and finally reaches the Buneman instability threshold (calculated for the initial, relatively 
low temperature T ej0 ~ 2 • 10 6 K). The subsequent development of the instability leads to 
anomalous turbulent heating (Sagdeev 1967; Biskamp & Chodura 1973; Bychenkov et al. 
1988; Kingsep 1991), which raises the electron temperature faster than the ion temperature. 
At some point, the electron temperature becomes much higher than the ion temperature 
and the ion-acoustic instability is excited. Such a transition from the Buneman regime to 
the ion-acoustic regime has in fact been studied previously (see Bychenkov et al. 1988 and 
references therein). In magnetic reconnection context, a possibility of this transition has 
recently been discussed briefly by Roussev et al. (2002). One complication, however, is that, 
since the plasma is constantly moving through the vicinity of the neutral point, then, in order 
to sustain the IAT, one needs to pre-heat the electrons (relative to the ions) in every fluid 
element that is just entering the diffusion region. This means that some sort of anomalous 
heat leakage across the magnetic field is probably needed. Whether and how this can be 
achieved is a difficult question, which falls outside of the scope of the present study. 
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5. Conclusions 

In this paper we have presented a model of magnetic reconnection in the presence of a 
current-driven enhanced anomalous resistivity. This is a very simplistic, crude model that 
aims at predicting the qualitative behavior of the system and the scaling of the reconnection 
rate with various plasma parameters, while treating numerical factors of order one only very 
approximately. 

In this model we have combined the following three ingredients. The first one is the ob- 
servation, derived from several recent resistive-MHD numerical simulations (Erkaev et al. 
2000, 2001; Biskamp & Schwarz 2001), that whenever the resistivity is strongly localized, 
the reconnecting system will develop a Petschek-like configuration, with the width of the 
inner diffusion region of the order of the resistivity localization scale. The second ingredient 
of our model is the Sweet-Parker model (Sweet 1958; Parker 1957, 1963) for the diffusion 
region of a Petschek configuration (Petschek 1964). Finally, the third ingredient is a physi- 
cally realistic model for a current-driven anomalous resistivity expressed as a function 
which exhibits two characteristic features. The first feature is a sudden jump of rj from a 
small collisional value i] to a much larger value rji as soon as j exceeds a threshold j c . The 
second feature is a subsequent linear growth rj oc j for j > j c . This choice is motivated by 
the theory of anomalous resistivity due to ion-acoustic turbulence, which has been developed 
in detail over the last 40 years (see, e.g., Bychenkov et al. 1988). 

Note that the anomalous resistivity function adopted in this paper becomes very sensitive to 
electric current density when the latter exceeds some threshold value j c ; this makes it possible 
for the resistivity to be enhanced only in a small region, which, in turn, leads to the devel- 
opment of a Petschek-like configuration. Thus, our model is characterized by a reconnection 
rate that is enhanced (with respect to the classical, collisional-resistivity Sweet-Parker rate) 
by a combined action of anomalous resistivity and of the Petschek mechanism. It is impor- 
tant to realize that the role of anomalous resistivity in the acceleration of the reconnection 
process is two-fold: in addition to its direct action (lowering the global Lundquist number 
S = VaL/tj), it accelerates reconnection indirectly, by turning on the Petschek mechanism. 

The width of the inner diffusion region of the Petschek model, and thus the resistivity 
localization scale, are determined self-consistently when all the ingredients of the model are 
taken into account. 

Based on our stability analysis of two possible Petschek-like states, we predict that the system 
will evolve towards a certain stable Petschek-like configuration. This stable configuration is 
characterized by the central current density jo an d the central resistivity r)(jo) exceeding j c 
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and 771, respectively, by a finite factor of order one. The reconnection velocity then scales as 
V Tec ~ Va/S* = T)i/5 C , where 5 C = cB /4:7ij c is the critical thickness of the layer. 

We then consider (in § 3) the case of anomalous resistivity due to ion-acoustic turbulence, 
as an important specific example. We derive very simple expressions for the parameters of 
the reconnection system (e.g., the width A of the diffusion region, reconnection velocity V Tec , 
and the reconnection time scale r rcc ) in terms of the basic plasma parameters n e , T e , Tj, 
and B . 

Finally, in § 4, we apply our model to the solar flare environment. We note that reconnection 
process will lead to significant electron heating, so that the electron pressure at the center of 
the reconnection layer may become comparable to the pressure of the reconnecting magnetic 
field outside the layer. Based on our model, we obtain typical reconnection times of order 
10 2 — 10 3 sec; this is short enough to explain the very fast time scale of impulsive flares. 
We note however that, as a result of the plasma heating inside the reconnection layer, the 
thickness 8 of the diffusion region quickly becomes comparable to, or even smaller than, the 
ion skin-depth, di = c/u pi . At these scales, new physical processes, described by Hall MHD, 
may come into play (Drake et al. 1994; Biskamp 1997; Bhattacharjee et al. 2001) even 
before the IAT develops and anomalous resistivity becomes important. 

I am grateful to S. Boldyrev, H. Li, R. Kulsrud, and R. Rosner for some very helpful discus- 
sions and interesting comments. I would like to acknowledge the support by the NSF grant 
NSF-PHY99-07949. 
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Fig. 1. — A schematic drawing of the inner diffusion region in the Petschek model. 
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Fig. 2. — The anomalous resistivity model adopted in this paper. 
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Table 1: Values of some basic plasma and reconnection region parameters in the typical 
coronal environment (column III) and in the flare environment (column IV). All values are 
calculated for pure hydrogen plasma (Z — 1, r/jj = m p ) and for ( = 2.0 (£ = e/( ~ 1.36); 
r rcc is calculated for L = 10 9 cm. 



Parameter 


Expression 


Typical Coronal Value 


Typical Flare Value 


T e ,K 


T 

-<■ e 


2- 10 6 


3. • 10 7 


Tii K 




2- 10 6 


3. • 10 6 


n e , cm -3 


n e 


1 • 10 9 


1 • 10 10 


-Bo, G 




100 


100 


f s , cm/ sec 


^T e /m p 


1.3 • 10 7 


5- 10 7 


Va, cm/ sec 


B /y/4irn e m p 


6.9 • 10 s 


2.2 • 10 8 


UJpg j SGC 


y/4irn e e 2 /m e 


1.8 • 10 9 


5.6 • 10 9 


(^Jpj, ; S6 C 




4.2 • 10 7 


1.3 • 10 8 


^~^g . SGC 


eB /m e c 


1.8 • 10 9 


1.8 • 10 9 


5 SGC 


ZeB§jm v c 


9.6 • 10 5 


9.6- 10 5 


Ai) e , cm 


A/T e /47m e e 2 


0.31 


0.38 


A£>j, cm 


y/Ti/Airriee 2 


0.31 


0.12 


d e , cm 




17 


5.3 


dj, cm 


c/uj pi 


720 


230 


j c , cgs-units 


2.14 en e t> s 


1.3 • 10 7 


5.1 • 10 8 


5 C , cm 


cB /47rj c 


1.8 • 10 4 


470 


771 , cm 2 /sec 


1.4 • IO- 3 (TJT^/Ujpe 


7.1 • 10 s 


2.2 • 10 9 




8im e T e /B$ 


7. • 10~ 4 


0.10 




2 • 10 4 (77T e )Vk/cV& 


1.7- 10 4 


45 


A, cm 




4- 10 7 


2.6 • 10 3 


r A (A), sec 




0.057 


1.2 • 10" 5 


Kec, cm/ sec 




1.6- 10 5 


1.9- 10 7 


E, cgs-units 


47rC 2 j c ^i/c 2 


5.2 • 10" 4 


6.4 • 10" 2 


7"rec 5 SeC 


L/V TCC 


6.3- 10 3 


52 



